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Abstract 

In this paper we study the evolution of multiple fluids with different constant densities 
in porous media. This physical scenario is known as the Muskat and the (multi-phase) 
Hclc-Shaw problems. In this context we prove that the fluids do not develop squirt 
singularities. 

1 Introduction 

We consider the dynamics of the two interphases between three incompressible and immiscible 
fluids in a porous medium without surface tension. The free boundaries are prescribed by the 
jump of densities between the fluids. Although the present paper is devoted to the evolution 
of two interphases in M 3 , the same approach can be performed when additional phases are 
considered in M. 2 or M 3 . The governing equations for the 2D incompressible porous media are 
identical with those modeling the dynamics in a Hele-Shaw cell (see |17j). 

The precise formulation of this problem is as follows pQ: the scalar density p = p(x,t) 
of the fluid is convected by the incompressible velocity flow u = {u\(x,t),U2{x,t),u^(x,t)) 
which satisfies Darcy's law i.e. 

pt + u ■ Vp = (Conservation of mass) 

V • u = (Incompressibility) (1) 

u = —Vp — (0, 0, p) (Darcy's law) 

where the scalar p = p(x, t) is the pressure and the acceleration due to gravity is taken equal 
to one to simplify the notation. 

Darcy's law yields the velocity in terms of the density by means of singular integral 
operators as follows: 

u(x,t) = PV ( K(x-y)p(y,t)dy - ^ (0,0, p(x)) , xGR 3 , (2) 

where the kernel K is given by 



K W = ~A 
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The integral operator is defined in the Fourier side by 

~tc\ £l+£iw t N 

U ^ = 1^|2 — )/=»(0 

that shows that the velocity and the density are at the same level in terms of regularity. 
The fluid is characterized by three different constant values of density p 1 < p 2 < p 3 

{p 1 in S7 1 = {x 3 > f(x 1 , x 2 ,t)}, 
p 2 in tt 2 = {f(xi,x 2 ,t) > x 3 > g(xi,x 2 ,t)}, 
p 3 in n 3 = {g(xi,x 2 ,t) > x 3 }, 

where f(xi,x 2 ,t) > g(xi,x 2 ,t). The moving surfaces 

Sf(xt,x 2 ,t) = {(xi,x 2 ,x 3 ) £R 3 : x 3 = f(xi,x 2 ,t)} 

S g (xi,x 2 ,t) = {(xi,x 2 ,x 3 ) £ R 3 : x 3 = g(x l7 x 2 ,t)} 

have the property (see below in section 2) that they can be parameterized as a graph for all 
time. Since the flow is incompressible the velocity of each interphase is continuous in the 
normal direction to the moving surface. Moreover, from the formulation it implies that the 
pressures are equal across the surfaces. The system is highly non-local in the sense that the 
equation for the moving surfaces involves singular integral operators and they are coupled 
together. Within the formulation we can recover the dynamics of a single interphase by 
taking p 1 = p 2 or p 2 = p 3 , which has been shown to be well-possed in the stable scenario 
with a maximum principle (see [7J and [8]). This case is known as the Muskat problem [16] 
(in 2D is also known as the two-phase Hele-Shaw) which has been broadly studied in [4], 
[TO] . [TH] . [7J, [8] and reference therein. 

The aim of this paper is first to show that for multiple interphases, in the stable case 
(p 1 < p 2 < p 3 ), the system is well possed in a chain of Sobolev spaces. In the unstable case 
(p 1 > p 2 or p 2 > p 3 ) there is Rayleigh- Taylor instability [7J and the system is ill-possed. 
Secondly we rule out a squirt singularity in the three phase system, i.e. that both interphases 
can not collapse in such a way that a positive volume of fluid between the interphases it gets 
ejected in finite time. Lets assume that both surfaces collapse at time T in a domain D such 
that 

lim [f(xi,x 2 ,t) — g(xi,x 2 ,t)] = for (x\ - x\) 2 + (x 2 - x 2 ) 2 < a 2 

where the constant a > and the point x = (x\,x 2 ,x 3 ) 6 D are fixed. Consider the domain 
Cl(t) denoted by 

tt(t) = {(xi,x 2 ,x 3 ) : (xi - xx) 2 + (x 2 - x 2 ) 2 < (R(t)) 2 , g(xi,x 2 ,t) < x 3 < f(xi,x 2 ,t)} 
with < R(t) < a and a section of its boundary given by 

S(t) = {(x 1 ,x 2 ,x 3 ) : (xi - xi) 2 + (x 2 - x 2 ) 2 = (R(t)) 2 , g(x 1 ,x 2 ,t) < x 3 < f(x 1 ,x 2 ,t)}. 
Then by using the divergence free vector field u it follows that 

— Vol n(t) = / [R'(t) -u-u]d (Area) 
dt Js(t) 
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where v is the unit normal to S(t). If the integral J \u\loocH is bounded then we can choose 
a time to £ [0, T) and take 

1 /- T 

.R(t) = - a — / | u| loo dr for to < t < T 

such that < < a for all t G [*o,^). Consequently ^ Vol 0(f) > for all t £ [t ,T) 
which prevents a collapse (squirt singularity) forming in between the moving surfaces (for 
more details see [5]). For a general n-dimensional definition of a squirt singularity see [6]. 
Therefore, the third part of the paper is devoted to show a bound of the velocity of the fluid in 
terms of C 1 ' 7 norms (0 < 7 < 1) of the free boundary. The estimate is based on the property 
that in the principal value ^) the mean of the kernels K are zero on hemispheres. This extra 
cancellation was used by Bertozzi and Constantin [2] for the vortex patch problem of the 2D 
Euler equation to prove no formation of singularities. For this system the convected vorticity 
takes constant values in disjoint domains and is related with the incompressible velocity by 
the Biot-Savart law (see |14| for more details). Let us point out that the system we are 
dealing with is a more singular one. Also, we quote the work [15] of Mateu, Orobitg and 
Verdera where this cancelation was used in quasiconformal mappings theory. 

Finally we would like to emphasize how the character of the kernels becomes crucial here 
since for analogous active scalar models we can not obtain this result. For the 2D surface 
quasi-geostrophic equation (SQG) [3] the kernels are odd (Riesz transforms) and for the 
patch problem [12] the velocity is not in L°°, it is in BMO (see [19] to get the definition and 
properties of the BMO space). Furthermore in the case of regular initial data for SQG and 
the system ([1]) the problem is also open [9]. 

The structure of the article is as follows. In section 2 we derived the equations of both 
interphases that are coupled together and in section 3 we show that the system is well-possed 
in Sobolev spaces. Finally, in section 4, we give a proof of boundedness of the velocity in 
terms of the smoothness of the interphases. 

2 The evolution equation of the moving surfaces 

The goal is to obtain the dynamics of a fluid that takes three different constant values of 
density p , p 2 and /? 3 as follows 

{p 1 in {x 3 > f(x,t)}, 
p 2 in {f(x,t)>x 3 >g(x,t)}, (3) 
p 3 in {g(x,t) > x 3 }, 

with f(x,t) > g(x,t) for all x = (xi,X2) £ M 2 . Hence 

Vp= (p 2 -p l ){d x J,d x J-l)8{xz -f(x,t)) + (p 3 -p 2 )(d xl g,d X2 g,-l)8(x 3 - g{x,t)), 

where 5 is the Dirac distribution defined by 

< h5,r] >= / h(x)r](x, f(x,t))dx, 
Jr 2 
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for rj(x%, X2, £3) a test function. 

For a divergence free velocity field, Darcy's law provides 

u = (d xl A~ 1 d X3 p,d 31iS A~ 1 d !B3 p, —d xl A~ 1 d Xl p - d X2 A~ 1 d X2 p), (4) 

and therefore 



u( Xl , X2 , X3 ,t)=-^PV f (I/I ,^V/(x -y,t).y) 



^ Jr* [\y\ 2 + (x 3 -f(x-y,t))^ 
p 2 - p 1 pv f (yi,y2,Vg(x -y,t) - y ) ^ 



(5) 



47r A 2 [|y| 2 + (^3 - g(x - y, t)) 2 f/ 2 

for f(x,t) / 23 / fi>(x,i). Taking 2:3 -> /(a?,i) yields 

u(x f(xt) t)= _ P 2 -P 1 pv f (yi,y2,Vf(x-y,t)-y) 

( ' /( ' j ' j 4vr FV L [\y\' + (f(x,t)-f(x-y,tm^ dy 

p 3 - p 2 m/ f (yi,y2,Vg(x -y,t) -y) 



4vr 7 K2 N 2 + (/(x,t)- ff (x-y,t))2]3/2 

without considering tangential terms. In any case the evolution of the surfaces are given by 
the normal velocity, the tangential terms are related with the parameterization of the free 
boundary [13]. Consequently 

(x,f(x,t)) t ■ (-d x J,-d x J,l) = u(x,f(x,t),t) ■ (-d x J,-d x J,l), 

and therefore 



ft t\= P 2 -P 1 py f (^f^^)-Vf(x-y,t))-y) 

P 3 ~ P 2 pv f Vf{x,t) -y - Vg(x -y,t) -y 
4vr J R2 [\y\ 2 + (f(x, t) - g(x - y, i)) 2 ] 3 / 2 V ' 



(6) 



In a similar way we obtain 



I A - P 3 -P 2 py f (Vg(x,t)-Vg(x-y,t))-y) 
9t[ , ) 4vr ™ UM 2 + {9<*,t)-9<*-V,t)) 2 ]> /a 
P 2 -P 1 pv f Vg(x,t)-y-Vf(x-y,t)-y 
4vr V [\v\ 2 + (<l(.x,t)-f(x-v,t))*]W V - 



(7) 



This coupled system describes the evolution of the moving boundaries Sf(x,t) and S g (x,t). 

3 Well-posedness in the stable scenario p 1 < p 2 < p 3 

Let us define the function d(f,g)(x,y,t) by the formula 

which measure the distance between the contours / and g. Therefore we consider f(x, t) 
approaching the value as \x\ — > 00 to avoid that the surfaces collapse at infinity. 
The section is devoted to prove the following theorem: 
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Theorem 3.1 Let f (x) - C tXll gQ{x) £ H k (R 2 ) for k > A, d(f ,g ) £ L°° and p 3 > p 2 > p\. 
Then there exists a time T > so that there is a unique solution to ([6]) and ([7|) given by 
(f(x,t),g(x,t)) where f(x,t) -C Ml g(x,t) E C 1 ([0, T];H k (M 2 )), f(x,0) = f (x) andg(x,0) = 

Proof: We shall show that the following estimate holds 

jE(t) < C(E(t) + If (8) 

for universal constants (C,p) and the function E(t) defined by 

E(t) = ||/ - C^Ws^t) + \\g\\ m (t) + \\d(f,g)\\ L oo(t). 

Applying standard energy estimates argument permit us to conclude local existence. This 
is based on introducing a regularized version of the system (|6|) and ([7]) that allows to take 
limits satisfying uniformly the a priori bound (jS}. 

To simplify the exposition we shall consider p 2 — p 1 = p 3 — p 2 = Air and k = A, being 
the rest of the cases analogous. Some of the terms can be estimated exactly as in [7] and 
therefore we shall show below how to deal with the different ones. 

We have 

5>-a*iii;C>=/ uv-cjpvf (v/w-v/(p,))., 

2 eft J R 2 J R 2 [\y\ 2 + (f(x) - f(x - y)) 2 ] 3 / 2 



+ / (f(x) - C oc )V/(x) • PV [ 

- / (f(x) - C^PV [ 
Jr 2 Jr- 



y 



[\y\ 2 + (f(x)-g(x-y)) 2 f/ 2 
Vff(x -y)-y 



\y\ 2 + (f(x)-g(x-y)) 2 f/ 2 
= h + h + h. 

For I\ we decompose further: I\ = J\ + J 2 + J3 where 

1 f f 1,1 \ r, , (Vf(x)-Vf(x-y))-y 



<i 



J2 = [ ( f(x) - C^PV [ 

J\v\>l J\y\>l 



V/(x) • y 



|v|>i ' ' J\y\>i[\y\ 2 + (f(x)-f{x-y)) 2 fl 2 



dydx 



and 



Since 



J 3= -/ (f(x)-c 00 )pv [ rrrZfv f(' V ^213/2 ^- 
J y \>i J\y\>i [\y\ 2 + (f(x) - f{x - y)YYi 2 

d x J(x) - d x J(x -y)= f Vd x J{x + (s - l)y) -yds, 

Jo 
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we get 

t <r f\ f I l-i f \f(x)-Coo\\V 2 f(x + ( S -l)y)\ 

Ji< ds \y\ / ,, 2 , 2W2 dxdy < C f - Coo L 2 Vf H i. 

Jo J\ y \<i Jm.* [1 + ((/(x) -f(x-y)) 2 \y\ 2 Y' 

In the term J 2 integration by parts yields 

2 J| y |>i Jr2 + ((/(x) - /(x - y)Yfl 2 

-2j\ y \J yl ) °°' ' [l + ({f(x)-f(x-yW\y\~^ ^ 

<C\\f- GoolliallV/Hioo <C\\f~ Coo|l! 2 ||V/||^. 

In a similar way for J3 we have 



Js= [ I (f(x)-C 00 )(f(x-y)-C ( 
J\v\>l JK 2 



2\y\ 2 -(f(x)-f(x-y)) 2 
\ y \>iW J ^ J ^ 0O,KJK " ^ C °'[\y\ 2 + (f(x)-f{x-y)Yf^ 



dxdy 



J\ y \>i Jri [\y\ 2 + (f(x) - f(x - y)) 2 fl 2 
- / / (/(*) - Coo)(f(x -y)- Coo) ^ ? f( rr^dxdaiy) 

J\ y \=i Jrz [\y\ 2 + (/w - f(x - y)rr 2 

^cdiv/iUoc + ijuz-Cooiii,. 

The term I2 is estimated as follows 



i*l I um - c^ {lix > - - »»W'> - »» • » 



2ju*jJ f{X> COO) [Is/I 2 + (/W - S (* - s)) 2 ] 5/2 

= 3 / / dxdy + 3 / / (ixdy 

/ 2 < § / |2/|~ 3 / |/0r) - CoopdV/^l + |V 5 (x - y)\)dxdy 

+ | / / l/^-CoopdV/^l + IV^x-y)!)^/^)^^)) 3 ^ 

< c||/ - Cooii^div/iUcc + ||v<7|| L oc)(K/, ^iiicc + 1) 

For I3 we write 

r f f tti \ r ^ t ^ 2|y| 2 - (/(x) - g (x - y)) 2 

*»=/ / (/(x)-Coo)(/(x-y)— — — — W^dxdy 

Jrz Jm.2 [\y\ 2 + (fix) - g(x - y)y\ b / 2 

+ f I (f(x)-C )o(x y) 3(/(x) " 9{X ~ y)){Vf{x) ~ V9{X ~ V)) • V - dxdy 

+ LL {m °° M V) l\y\ 2 + (m - g(x - y)) 2 W 2 V 
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and proceeding as before we find 

h < C\\f - OolMlslMllV/IU- + \\Vg\\ L °o + l)(M(/ j5 )||ioo + 1). 
Thus Sobolev inequalities yields 

j t \\f-C O0 \\ 2 L2 (t)<C(E(t) + lf + \ 
The analogous estimate for g is 

±\\ g f L2 (t)<c(E(t) + iy + \ 

To estimate the higher order derivative, we consider the quantity 

^|ll^ 4 1 /ll! 2 (t)=/4 + / 5 +/6, 



where 



h= d f(x)d [PV — — ^T^dyjdx, 

Jr 2 v Vr 2 [\yr + U\ x ) - f( x - y)) \ ' ' 

h = L a -- /( ^ ( PV L wHm-l(*-vm^ > x 



and 



Jr 2 v Jr 2 [\y\ 2 + (/0M) -g(x-y, t)) 2 ] 6 ^ 

The estimates for I± are obtained in [TJ. 

We split I5 and consider the most singular terms which are 



J*= I dij{x)Vd A x J{x).PV I 

JR2 JRi 



dydx, 



[\y\ 2 + (f(x)-g(x-y))*]W 



and 



„ , .,, ,, , 1 ^f{x)-y{f{x)-g(x-y)){&^J{x)-d^ 1 g(x-y)) 



il 2 il 2 ' 



[\y\ 2 + (f(x)-g(x-y)) 2 } 11/2 



For integral J4 we find that 



Ji ~ 2 L L l9xJ[x)l [\y\ 2 + (fix) - g(x - y)) 2 f' 2 dxdy 
< CH^/H^CIlV/IUoo + \\Vg\\ L ~)(\\d(f,9)\\l~ + 1). 
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For J5 it follows 



■h 



\v\ 



+ 



|y|<i 



\d x J(x)\\Vf{x)\{\d x J{x)\ + \d xi g(x - y)\)dxdy 



\dij{x)\\Vf{x)\{\dtj{x)\ + \dt ig {x-y)\)(d{f,g)(x,y)fdxdy 



< CllC/II^IIV/ll^dl^/H^ + \\dt ig \\ L 2)(\\d(f,g)\\U + I)- 
Similarly we obtain 



J fi = C 



\y\- 



bl>i 



+ C 



\y\<lJR2 



K f(x)\\Vf(x)\(\d x J(x)\ + \d Xl g(x - y)\fdxdy 



\d x J(x)\\Vf(x)\(\d Xl f(x)\ + \d Xl g(x-y)\)Hd(f,g)(x,y)fdxdy 



< cii^/iu.iiv/ii^di^/iiioc + H^IU-nK/^llioc + 1) 



and finally 



Jt \\d x j\\l,{t)<c{E{t) + iy+\ 



Analogously we get 



j t \\d$ 2 f\\Ut)<c(E(t) + iy+ 1 , 



and therefore using ([U]) it follows 

d 



dt 



||/ -CUff^t) + 1) 



p+1 



In order to estimate the H 4 norm of g we proceed as for the surface / to obtain 

±\\gf m (t)<C{E(t) + lf+\ 



(11) 



(12) 



Next, we analyze the evolution of the distance between the surfaces. For the quantity 
d(f,g)(x,y) we have 



^(/, g)(x,y,t) 



(fix) - g(x - y, t))(f t (x) - g t (x - y)) 
[|y|2 + (/(x)-<7(x-y))2]3/2 



(13) 



< (d(f,g)(x,y,t)f 



till/ 



>(t) + \\g t \\L~(t)). 



We now estimate ||/t||i«>(t), being equivalent to control ||<?t||i,°°(i). We consider for f t (x,t) 
the following splitting: 



I 7 = PV 



(yHx,t)-Vf(x-y,t))-y 
[\y\ 2 + (f(x,t)-f(x-y,t))^ 

8 



dy, 



r = PV [ Yli^tl 
8 J*> [\y\ 2 + (f(x,t)-g(x-y,t))*]W V > 

T = _ PV [ Vg(^ -y,t)-y , 

L [\y\" + (f(x,t)-g(x-y,m^ ^ 



Let us decompose further: I-j = J-j + Jg + Jg where 



T = f (Vf(x,t)-Vf(x-y,t)).y 
7 J lyl<1 [\y\* + (f(x,t)-f(x-y,t)r]m y ' 



J 8 = Vf(x,t)-PV I 

J\v\ 



y ,1 

ay 



and 



Hence 



l\y\>l [\y\ 2 + (f(x,t)-f(x-y,tm^ 

J9 = -pv f — — " ~-^<iu- 



V/(x — y,t) ■ y 
',|>i [\y\i + (f(x,t)-f(x- y ,t)r]^ 



\Jr\<C [ lyr^yllVVlU- < C\\V 2 fl 

J\v\<l 



'\y\<i 
We rewrite Jg as 

T=Vf(rt) f V ^-^ + (f(^)-f(x-y,t)f\yn^ 
8 n ' j '4|>iH 3 [l + (f(x,t)-f(x-y,t)r\y\-^ dV > 

and considering the function Q(a) = [1 + a 2 ] 3 / 2 , the mean value theorem gives 

T = XIf( t\ f JL -3[l+/3 2 ] 1/2 /3|/(x)-/(x-j/)| 

8 V/l ' ' ' 4|>i M 4 [1 + (f(x,t) - f(x - y^M- 2 ? 12 ^ 

where < (5 < \ f(x) - f(x - y)\\y\~ l . Therefore 

|J" 8 | < C||V/|Uoo||/ - CooHioo. 

In a similar manner as for J3, integrations by parts in Jg yields 

\M ^cdiv/iUoo + ijuz-CooiUco. 

For Is we split 

h = PV ( dy+ I dy = K l +K 2 . 

J\y\>l J\y\<l 

We deal with K\ as we have done with Jg to get 

\K X \ < ||V/|U=o(||/ - CoolUco + Coo + HflrlUoo). 
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For K2 it is easy to obtain 

\K 2 \<\\Vf\\ L ~\\d{f,g)\\l»°. 

We control Ig as the term ^3 and therefore 

|/ 9 | < CUffHiocdlV/ll^ + \\Vg\\ L - + l)(||d(/,<7)||i« + !)■ 
The above estimates together with (fT3j) yields 

^d(f,g)(x,y,t) < d(f,g)(x,y,t)C(E(t) + 1) P , 

and integrating in time 

d(f, g)(x,y,t + h)< d(f, g) (x, y, t) exp ( C(E(s) + lfds 
for h > 0. Hence 

IK/, 0) IU- (* + /i) < IK/, <?) |U=o (t) exp ( jf ' + " C(E( S ) + lfds 
The above estimate applied to the following limit: 

d„. r , x„ m l|d(/,g)l|L°°(t + ^)-||d(/,g)IU°°(t) 

5 l|d(/^)||L-(t) = 

allows us to get finally 

j t \\d(f,g)\\ L ~(t)<C(E(t) + lT. (14) 



Combining estimates (jlip . (|12p and (|14p we conclude that (|8|) holds for universal constants 
C and p. 

For uniqueness we proceed as in [7] which leads to the desired result. 

4 Bound for the fluid velocity: even kernels 

In this section we prove the following lemma: 

Lemma 4.1 Let f,g be solutions of the contour system (J6|) and ([7]). Then the velocity of 
the fluid satisfies the following bound 

NU- < C(l + - + - ln(l + || V/Hcr + ||V 5 ||c-r) 

V T T (15) 
+ ln(l + ||/ - Coolie + Coo + ||V/|| L2 + llpHi- + ||V 5 || L2 )) , 

where < 7 < 1 and the constant C = C(p 1 ,p 2 ,p 3 ) depends on p , p 2 and p 3 . 

10 



Proof: We shall denote x G M 2 and x = (x,x 3 ) G M 3 . In order to acquire the above 
inequality, we can split p as follows: 

P = p 1 Xn+(t) + P 1 x^(t) + P 2 xn2(t) + (p 3 ~ P 2 )xns(t) + P 3 Xn-(t) 

where 

n + (t) ={x 3 > M}, ST(t) = {x 3 < -M}, 
n\t) = {M>x 3 > f(x,t)}, n 2 (t) = {f(x,t) >x 3 > —M} 

and 

n 3 {t) = {g(x,t) >x 3 > —M} 

for M = ||/ - Coo Ik- (t) + Coo + \\g\\L°°(t) + 1. 
From Q we check that 

(d X3 A~ 1 d Xl ,d X3 A~ 1 d X2 ,-(d Xl A~ 1 d Xl + d X2 A~ 1 d X2 ))(xn±(t)) = 0. 
Thus the Fourier transform yields 

u = T(p 1 xn^t) + p 2 Xn2(t) + (p 3 ~ P 2 )xnz(t)) 

2 ( 16) 

- ■^(0,0,P 1 Xtt 1 (t) + P 2 XnHt) + (P 3 ~ P 2 )xn?{t)), 

where 

i /■ 

7TT r- Tn>3 



T(h)(x) = ^PV [ K(x - y)h(v)dy, x G 



and 

K{x) = (&^ 2x "-^- xl 
\ \x\ b \x\° 

Then we have to deal with 

p 1 T(xnm)), P 2T (Xn2(t)) and (p 3 - p 2 )T{x&(t))- 

We will consider T(xn 1 (t)) since the other terms are analogous. We proceed as in [2]. The 
key point is that the kernels in T are even. Therefore in the principal value the mean 
of the kernels are zero on a hemisphere. We consider the three different coordinates of 
r(xoi(t)) = (Ti,T2,T 3 ). Then for a fixed x G M 3 we have 

4vr J n i(t) \x-y\ b 

We take a distance 5 given by 

1 

8 



Zil + WVfWc^ + WVgWc-,) 1 ^' 
Then if d(x,n l (t)) >Swe consider O x (t) = U l (t) U ?7 2 (t) for 

U\t) = n 1 n {(xi - yi) 2 + (x 2 - y 2 ) 2 < L 2 }, 
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and 

u\t) = n 1 n {{ Xl - yi f + (x 2 - y2 f > l 2 }, 

where L = 2(1 + ||/ - C^U*. + + \\Vf\\ L i + \\g\\ L oo + \\Vg\\ L 2). 
The splitting T x = h + 7 2 for 



3 f (»,-»,)(*,-»,) 

3 dt/ /" (^i -z/i)(^3 -ys) 



47T Jtj: 



lu 2 (t) \x-y\ b 
gives 

\h{x)\ < C / r _1 dr < Cln(\/2L/<5) 
Js 

We write I 2 = lim^+oo such that 



dy 



where 



4^ Ju\(t) \x - y\ 6 



U 2 R (t) = ft 1 n {L 2 < ( Xl - yi f + (x 2 - y 2 ? < R 2 }. 
Then integration by parts gives 



t r 3 f (x 3 - yi) 
h = -7-1 1= — =Tn-nido{y), 
4 7T Jdul(t) \x-y\ 6 



ldu 2 R (t) \x-y\ z 

and therefore 



if 



-! 

4vr J L< \ X . 



d Xl f(y)(x 3 -f(y))dy 3_ [ 2 * f M (x 3 - y 3 )Ly' 1 d6dy 3 

1 - 'l<\ X - V \<R (\X ~ y\ 2 + (X 3 " f(y)) 2 ?' 2 4VT J Q J f{x+Lyl) (L 2 + (x 3 - V ,) 2 f/ 2 



3 f 2n f M (x 3 -y 3 )Ry' 1 d0dy 3 



4vr7 Jf( x+ Ry>)(R 2 + (x3-y3) 2 ) 3 / 2 
for y' = (cos 6, sin 0). By the Cauchy-Schwarz inequality we find easily that 

lT R {< 3 \\d Xl f\\ L , 3M , 3M 

We take now R to infinity to get \I 2 \ < C for C & universal constant. 
If d(x, ^(t)) < <5 we set 7\ = Ji + J 2 where 



4vr 7ni(t)nB,(x) k — 2/T 



3 DT7 /" (a?i -yi)(^3-2/3) 

and 



m = iipy [ 

4tt Jq 
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For Ji one can proceed as in [2\. It is clear that Jt(x) = if d(x, dti 1 ^)) = d>5. Therefore 
we take d < 5. We show the argument for the boundary of Q, l (t) given by x 3 — f(x,t) = 0. 
The planar section trivializes. Then we define the following set of directions 

S r (x) ={z<ER 3 : \z\ = l,x + rze n\t), d < r < 5} 

Also, we pick the point a £ M 3 on the boundary, a 3 — f(a, t) = 0, such that \a — x\ = d. We 
define the following hemisphere: 

S+(x) = {z G R 3 : \z\ = 1, (-V/(a,i),l) • z> 0}. 

For the directions z in the symmetric difference 

D r (x) = (S r {x) \ S + (x)) U (S + (x) \ S r (x)), 

we define the angles <p(z) by 

(-V/(a,i),l)-Z 



sin tp(z) 



y/l + \Vf(a,t)\2\z\ 



It is easy to check that simp(z) > gives x 3 + rz 3 — f(x + rz) < and smtp(z) < gives 
X3 + rz 3 — f(x + rz) > 0. Therefore 



simp(z)\ 



(-V/(a, t), 1) • (a - x) + _(-V/(a, t), 1) • (x + rz - a) 



Vl + |V/(a,t)|2 



Vl + |V/(a,i)| 2 : 



r 



(-V/(a,i),l) • (sc + rz-a) (a 3 - /(a,*)) - (x 3 + rz 3 - f(x + rz,t)) 



Vl + |V/(a,t)| 2 r 



+ 



v/l + |V/(a,i)| 2 ' 



< d ] |/(x + rz,t)-/(a,t)-V/(a)-(x+rz-o)| 
— r r 

and finally 

I sin <p(z) I < - + H V /HH^ + rz - a r 7 < rf + l|V/IM^+0 1+7 

y» 

< d + o-v llV/llc^^ + r^) < d + 2 ^i+7 +r i+7 ^ (1 + 27) d + /2r > s - 

— r r — r fty r - r 



5 ) 



The above estimate and the fact that 



>S+(x) 

yield for J\ the following bound 

-5 /•27T 



(x l - yi)(x 3 - y 3 ) 

— — da{y) = 

\x - y\ 5 



\Ji(x)\<— / / X{ v (z):zeD r m\ cos <p\dipd9- < — (l + 2T(l + l/ 7 )). 
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We write J2 = K\ + if 2 where 

(xi - yi)(x 3 - y 3 ) 



Kl {x) = £-PV f 



L (*)nB|(x) k-2/l 5 

and 



(xi - yi){x 3 - y 3 ) 



dy, 



dy. 



lu 2 (t)nB c s (x) \x-y\ 5 

Then one can deal as for I\ and I2 respectively to obtain the same estimates. 
We write the second coordinate as follows 



4vr J n i {t) {x-y^ 



to get an analogous bounds. 

For the third coordinate we have 



Ux) = ±pv [ ^-yrf-^-yf-^-y*) 2 ^ 

4vr J n i (t) \x-y\ 5 



and we can proceed as before. But for the term we have 



T R 



J_ f d ( {X3 ~ y3) )dy- 1 f f M d -M*) - 

1 f (n^)-f( T? ))A. 



4tt J L <\x- y \<R \x-y\ \x-y\ \x-y 

where P(a) = a/(l + a 2 ) 3 / 2 . The mean value theorem yields 

t r_ 1 /" nlf ^(f(y)-M)dy 



- 1 1 

47T ,/l,<|x-i/l<.R 



1 13 
'L<|x-?/|<R l x ~~ 2/1 

and using that \P'(a)\ < 2 we find easily < 2M/L < 1. 
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